Spatial Coherence of a polariton condensate in ID acoustic lattice 
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Several mechanisms are discussed which could determine the spatial coherence of a polariton 
condensate confined to a one dimensional wire. The mechanisms considered are polariton-polariton 
interactions, disorder scattering and non-equilibrium occupation of finite momentum modes. For 
each case, the shape of the resulting spatial coherence function g^ 1 ' (x) is analysed. The results 
are compared with the experimental data on a polariton condensate in an acoustic lattice from [E. 
A. Cerda-Mendez et al, Phys. Rev. Lett. 105, 116402 (2010)]. It is concluded that the shape of 
g' 1 ^ (x) can only be explained by non-equilibrium effects, and that ~ 10 modes are occupied in the 
experimental system. 



I. INTRODUCTION 

The fabrication of planar (2D) microcavities has pro- 
vided experimental access to the strong coupling regime 
between light and matter in a spatially extended system 
The excitations of the coupled modes are polaritons, 
quasi-particles which are bosons with a very small effec- 
tive mass, typically 10 -4 times that of an electron. This 
makes it possible to study quantum effects in a solid state 
system at a relatively high temperatures The short 
life time of the polaritons, ~ 5 ps means that the 
system has to be constantly pumped to observe Bosc- 
Einstein condensation (BEC). It can be achieved either 
resonantly [3] or non-resonantly 0| producing conden- 
sates with long coherence times [5|] and long range spatial 
coherence [3, 6j,[l2j]. These have been shown to exhibit in- 
teresting nonlinear and many-body phenomena, such as 
vortex formation [9|, llfj ] , solitons [8j and superfluid-like 
flow 0. 

In recent experiments (l2l. f]~3j| . microcavity polaritons 
have been confined in wire like geometries, to study the 
peculiar interaction effects of bosons in a ID system as 
well as a step to the realisation of polariton circuits 
flll |. The extra confinement was provided in two dif- 
ferent ways: b y e tching of the whole planar structure 
to form a wire [12J and by excitation of surface acoustic 
waves (SAW) to form a dynamic ID lattice [13]. It is 
the spatial coherence measurements in the latter setup 
that are the subject of the present paper. In the ex- 
periment, the same sample area was studied with and 
without the SAW, so it was possible to measure how the 
coherence length changed when the ID confinement was 
introduced. Without confinement, the entire 2D emis- 
sion spot was coherent, that is, the coherence length was 
the same as the spot size (l3l . With confinement it 
was observed that the coherence length in the direction 
perpendicular to the SAW wavefronts was reduced to the 
SAW wavelength, as the condensates confined in the in- 
dividual minima became decoupled. However, there was 
also a significant reduction in the coherence length mea- 
sured parallel to the wavefronts 13|. along the 'wires', 
which is less easily explained. 

In this paper we analyse possible mechanisms for the 



reduction of the spatial coherence along the ID wires. 
We consider the effect [l5j of polariton-polariton interac- 
tions, disorder, and finite momentum modes occupation 
in a harmonic trap on the first order correlation func- 
tion </W (x). The theoretical predictions are compared 
with the experimentally measured g^ (x) [lq | which has 
a Gaussian shape up to the noise floor. This shape can 
only be reproduced by the non-equilibrium occupation of 
a set of the finite momentum modes. Performing a more 
quantitative fit, using the occupation numbers as fitting 
parameters, we find that the number of the excited modes 
in the experiment [13] is ~ 10. 

The rest of the paper is organised as follows. In Section 
2 we formulate the model of the ID polariton system with 
different types of interactions. In Section 3 we calculate 
the spatial dependence of the coherence function (x) 
when finite momentum modes are occupied in a harmonic 
trap. Section 4 contains discussion of the Anderson local- 
isation due to disorder in a ID system. In Section 5 we 
discuss the effect of the short range polariton-polarion in- 
teraction on (x). Section 6 contains the fitting of the 
experimental data with the theoretical predictions from 
the previous sections. 



II. MODEL 

We consider one dimensional polaritons formed by the 
equal mixture of exciton and photon at the zero momen- 
tum following the experimental setup of (l3| . 

The motion of photons in z— direction is restricted by 
a planar microcavity which is formed by a pair of Bragg 
mirrors giving a small mass to the 2D photons Ck. Exci- 
tons are bound electron-hole pairs, they are bosons X^, 
which are confined to the same 2D plane by a quantum 
well, see the scheme in Fig. 1. The eigenmodes of a 
strongly interacting X^ and Ck are two hybrid particles, 
Lower and Upper polarions (LP and UP). We will con- 
sider only the LP branch as bosons with a parabolic 
dispersion and mass m^p assuming that the photons and 
the excitons are at resonance at k = 0, e. g. a LP around 
the zero momentum consists of equal mixture of the ex- 
citon and the photon Ok = {X^ + Ck) /v2- 

Bosonic nature of these polaritons allows a macroscopic 



occupation of the single zero momentum state that leads 
to a long-range phase coherence in the system. Access 
to this regime is not straightforward as the short lived 
polaritons do not have time to relax at low momenta due 
to a slow emission of the low energy phonons forming 
a relaxation bottleneck and forbidding formation of the 
truly equilibrium condensate. This difficulty can be over- 
come in a nonequilibrium regime under a strong pumping 
when the k = state is populated directly or indirectly 
by employing an optical parametric oscillator (OPO). 

In the OPO scheme the system is excited by the pump 
at a finite momentum k p and energy E p . Due to the 
polariton-polariton interactions the k p excitations can 
scatter fast to a low momentum. Above a threshold den- 
sity of the directly pumped polaritons the interaction be- 
comes strong enough to reach the zero momentum state 
k s = 0, signal, which also forms a significant population 
in the idler state kj due to the conservation of the en- 
ergy Ei + E s = E p and the momentum k, + k s = k p 
by the scattering process. Above the threshold the sys- 
tems locks only to these three states Q with a significant 
population of the signal state at k s = 0. This approach 
is more advantageous than the direct pumping as in the 
OPO regime the pump at a finite k p does not obstruct 
the emission from the condensate at k s = 0. 

The 2D polaritons are constrained further to a set of 
ID wires by an in-plane surface acoustic wave propagat- 
ing along y-direction, see the scheme in Fig. 1. The 
y-component of the momentum belongs to the lowest en- 
ergy subband leaving the only motion along a;-axis un- 
constrained. The polaritons are detected though the pho- 
tonic part which constantly escapes the system forming a 
spatial distribution of the electric field outside of the cav- 
ity that contains information about the polaritons inside 
the structure, 
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and can be registered by a detector. Here ipk (x) are the 
eigenmodes of the external potential in cc-direction, e. g. 
<fk (x) = e tkx if the there is no extra potential along re- 
axis, e/c is the dispersion of the electro-magnetic modes, 
and L is the length of the wire. 

We consider three possible mechanisms that can af- 
fect spatial coherence of free polaritons in the signal, a\ 
around fc s = 0. First is a confining potential originating 
from a finite size of the excitation spot. We assume that 
the single particle potential is harmonic, 
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Figure 1: Scheme of the microcavity (DBRs) with a quantum 
well inside (QW), and the surface acoustic wave (red elon- 
gated ellipses are the polariton ID wires). 



where F (y) is a function that vanished when \x — x'\ 3> I, 
i. e. F (0) = U and F (y > 1) = 0. And the third 
mechanism is the pair polariton-polariton interaction. It 
has a short range interaction potential as the polaritons 
have zero electric charge with the interaction strength Vb . 

We neglect the interaction of polariton with phonons 
mechanism for decoherence as it is ineffective at low mo- 
menta. The polariton-phonon scattering time is ~ 100 
ps [17j which is much longer than the polariton life time 
~ 5 ps. The latter is due to detuning chosen in the ex- 
periment [l3[ where the polaritons at k ss consist of 
the equal mixture of the exciton and the photon. 

To analyse the spatial coherence that is observed via 
the photonic part of the LP we consider the first order 
correlation function, 
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where N [. . .] is the normal ordering operator that elim- 
inates the vacuum fluctuations of the quantised electro- 
magnetic field and (. . .) is the average with respect to a 
state of the system. 

The eigenstates of the free polariton Hamiltonian are 
plain waves. The expectation value over the zero momen- 
tum state occupied by many polaritons gives an infinite 
range coherence, <?W (x) = 1. In the next three sections 
we analyse each of the possible mechanisms separately to 
see how they restrict this behaviour. 



III. FINITE MOMENTUM MODES IN A TRAP 



where the frequency u> defines the system size as L = 
\j2hj (rriLpuj). Second is a disorder potential that has a 
finite correlation length Iq, 
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Here we obtain the single particle eigenstates of polari- 
tons in the harmonic trap. Then we evaluate g^- 1 ' (x) as 
the expectation value with respect to a density matrix 
assuming a non-equilibrium distribution of the polariton 
occupation numbers around the zero momentum state. 
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The single particle Hamiltonian for a polariton in a 
harmonic potential, 



Substituting Eqs. ((HJ in Eq. (Q} and using the eigen- 
function Eq. ^ we obtain the first order coherence func- 
tion as sum over a set of finite momentum modes, 
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is the one of the quantum harmonic oscillator. It is a well 
known and solved model. 

Following the standard approach to the diagonalisation 
problem we search for a solution in form of a polynomial 
|19| . The eigenfunction problem for the model Eq. ([5]) 
leads to Hermite's differential equation which is solved 
by 
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where Hk(z) = (—l) k e z d k e z /dz k are the Hermite 
polynomials, L — y/2H/ (m^pco) is the size of the system, 
and k is the number of the energy level that describes a 
finite momentum state of the ID polariton. The comple- 
mentary eigenvalue problem can be solved via a Taylor 
expansion. As a result the eigenenergies of the finite mo- 
mentum states k ^ have a linear spectrum 



E k = uj I k + I 



(7) 



A many polariton state, which is created in a pumping 
process, is described by a density matrix. The expecta- 
tion value in the correlation function Eq. (0| has to be 
evaluated as an ensemble average, 
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where p is a density matrix. 

In a continuous wave experiment the signal is accumu- 
lated on the time scale of seconds. The phase coherence 
of a polariton condensate is limited by few hundreds of 
pico-seconds jj| which partitions the long accumulation 
interval into a large ensemble of many very short mea- 
surements. Thus the total signal is averaged over many 
realisations with different distribution of the polariton 
occupation numbers and the total number of polaritons 
which also changes slightly in time due to the power fluc- 
tuation in the pumping laser jHJ. As the phase coherence 
is lost between different realisations the density matrix, 
which describes such an ensemble, is diagonal in the rep- 
resentation of the polariton occupation numbers, 
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Here \noni . . .) is a state with no polaritons at the mo- 
mentum k = 0, tlx polaritons at the momentum k = 1, 
etc and Pn oni ... are the diagonal matrix elements that 
describe a non-equilibrium distribution of the boson oc- 
cupation numbers. 



(x) =22^k (x)ipk (- 



X)Pk, 



(10) 



k=0 



where p k = En rH... n kPn oni .../ E„oni. : . /Wi- are 

weight factors, e. g. in thermal equilibrium they are 
given by the Bose-Einstein distribution function. Note 
that pk are not completely equivalent to the occupation 
number of the polaritons, rik, due to the normalisation 
of the correlation function in Eq. (j4)) . 

When only a few low momenta modes are excited with 
probabilities that are approximately equal, pk ~ pi, the 
shape of the correlation function is almost Gaussian with 
a reduced system size L < L. gV-> (x) ~ e~ x l L . When 
many modes are excited with an arbitrary p^ the shape is 
arbitrary as polynomials in Eq. (JSJ) of a high order k 1 
can approximate a wide range of different functions. 



IV. ANDERSON LOCALISATION 

Polaritons are subject to a random potential which 
originates from two sources. One is disorder of the quan- 
tum well that influences the exciton component of polari- 
tons. And another is imperfections of the cavity mirrors 
that influences the photon component [U [22| ■ Here we 
consider a simplified model of the disorder characterised 
by a single correlation length Eq. (J3J) and we assume that 
the wave length of the surface acoustic wave is smaller 
then the correlation length of disorder potential in the 
2D heterostructure to model the strictly ID system. If 
the latter is not true the polariton wire becomes a quasi- 
1D system which would give a longer localisation length. 
In this section we also neglect interplay between disorder 
and polariton-polariton interactions which can lead to a 
many body metal-insulator transition above a finite tem- 
perature [23l Izij and consider the low temperature case 
only. 

Scaling analysis of the potential disorder problem (20| 
shows that conductance of a ID system is zero at low 
temperatures as all of the zero momentum states become 
localised when the system size increases to infinity. The 
lowest energy eigenstates of the single particle model, 
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are localised with the exponential tails \ip (x)\ ~ e(' x ''&, 
for example see a review |25j. Therefore the first order 
correlation function in Eq. (U) evaluated with respect to 
these states gives a finite coherence length with the same 
exponential tail, 



(12) 
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where, in general, the localisation length £ is not equal 
to the correlation length of the random potential l Q . 

Here we do not perform a more detailed study of the 
features that are specific to the disorder in the polariton 
systems, e.g. [2(|, but only use the exponential tails of 
the correlation function as a characteristic feature of the 
localisation mechanism. 



V. POLARITON-POLARITON INTERACTIONS 

The Mermin- Wagner theorem [18| states that a long 
range order in a ID system is absent due to any finite- 
range exchange interaction which results in a large long- 
range quantum fluctuations. In this section we discuss 
the manifestation of this general statement in the specific 
system of interacting ID polaritons. 

The exciton component of a polariton has zero charge 
but a finite dipole moment. Neglecting the photonic non- 
linearity, interaction between two polaritons has dipole- 
dipole nature thus it is short range. This interaction 
limits temporal coherence of the polariton condensate 
Q. The spatial coherence, within the model of density- 
density interaction with the delta-function interaction, 
was was analysed in details in (27l HI| • 

It was shown that the coherence length is finite which 
is manifested by the exponential tails of the coherence 
function, 

S « (re) ~ e -*/S (13) 

in accordance with the general Mermin- Wagner theo- 
rem [l8| that forbids an infinite range order in the ID 
systems. It was also shown in [28] that the coherence 
length decreases with increase of the interaction constant, 
l c ~ 1/Vb- For a typical GaAlAs microcavity the estimate 
of l c was given as few hundred micrometers in an optimal 
regime. 

Note that the functional form of the coherence function 
due to the polariton interactions Eq. (|13[) and due to 
disorder Eq. (j!2p coincides in the tails. 




Figure 2: Coherence function g ( ' measured along ID wire on 
the logarithmic scale. The thin dashed line marks the noise 
floor and the thick dashed line is an exponential function. 




Figure 3: Coherence function g ( ' measured along ID wire on 
the normal scale. Thick line is the numerical fit using the Eq. 
Q with pfe as the fitting parameters. Thin dashed line is the 
plot of Eq. Q with pt given by full triangles in Fig. [4] 



VI. FITTING OF EXPERIMENTAL DATA 

Here we analyse experimental data [l^, [l6| on the co- 
herence of polariton condensate along ID wires. 

In this experiment the 2D spot of the pumping laser de- 
fines the system size and, therefore, limits the coherence 
length of the polaritons. Then the microwave radiation is 
applied to confine the polaritons to a wire-like geometry. 
It was observed that the coherence length in the direction 
perpendicular to the wavefront was reduced from the sys- 
tem size to the wave length of the surface acoustic wave 
but, unlike it was expected, the coherence length in the 
unconstrained direction, parallel to the wavefront, was 
also reduced. 

We start from the qualitative analysis of the shape 
of the measured coherence function, g^, to identify the 



main mechanism that reduces the coherence length along 
the wires. The data on is presented in Fig. [5] on a 
logarithmic scale. The thick dashed lines marks an expo- 
nential function, which would correspond to the Ander- 
son localisation or polariton-polariton interaction mech- 
anisms from Sections 4 and 5, and serve as guides for 
the eye. Investigating both tails simultaneously, starting 
from the noise flow marked by the thin dashed line in Fig. 
[21 we find that the measured function has a super-linear 
character away from the maximum point [29]. Thus we 
conclude that the main mechanism is the non-equilibrium 
occupation of a set of the finite momentum modes from 
Section 3. 

Having identified the mechanism we perform a more 
quantitative fitting. We extract the 2D spot size from 
the data on </W which was measured without the surface 
acoustic wave [16]. Gaussian fit gives the system size 
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Figure 4: Occupation numbers p^ in Eq. (|10p to fit the exper- 
imentally measured g"' : full ellipsis - result of a numerical 
fit (think solid in Fig. [3}, full triangle - a set of p^ chosen by 
hand to occupy approximately 10 modes (thin dashed line in 
Fig. E). 



L = 16.3 fim. Then we fit the data on g'- 1 ' , which was 
measured with the surface acoustic wave, by the corre- 
lation function from Eq. f|10|) using pk as the fitting pa- 
rameters. The gradient descent method gives the thick 
full line in Fig. [3J The set of pk for the thick line, which 
is presented by filled ellipses in Fig. 21 shows that ap- 
proximately ten modes are occupied. 

This fitting of a single function with many parameters 
is not unique. Another set of pk that is chosen ad hoc to 
occupy approximately the same number of modes, trian- 
gles in Fig. 21 also fits satisfactory the experimental data, 
thin dashed line in Fig. [3] It is not possible to extract the 
set of pk from the data on g^ uniquely but characteristic 
number of the occupied modes is determinable. 

We also analyse distribution of p^ at thermal equilib- 
rium. Bose-Einstcin distribution, 
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where is the single particle spectrum of a harmonic os- 
cillator from Eq. Q and f3 is an inverse temperature, at 
low temperatures gives the distribution of the occupation 
numbers as p^ ~ 1/fc. It does not fit a satisfactory the 
experimental data. Thus we can conclude that polaritons 
are not in thermal equilibrium. This reflect the strong 
out-of-equilibrium nature of the constantly pumped po- 
lariton system. 



VII. CONCLUSIONS 

We have analysed different mechanisms that can limit 
the spatial coherence of a polaritonic condensate in ID 
wires formed by the surface acoustic wave and have com- 
pared their predictions with the shape of the first or- 
der coherence function g^ 1 ' (x) measured in [Ijj. From 
a qualitative analysis of the experimental data we have 
found that the main mechanism is the non-equilibirum 
occupation of a set of the finite momentum modes. Out 
of the three effects that we considered: Anderson local- 
isation and the polariton-polariton interaction give an 
exponential tails of g^ (x), the finite momentum modes 
occupation gives a Gaussian shape. In the experimental 
data the shape is Gaussian, up to the noise floor. 
Performing a more quantitative fit, using the occupation 
number as the fitting parameters, we have found that 
the number of the excited modes in the experiment (l3| 
is - 10. 

From the analysis in the present paper it is advised 
to use a static method of the polariton confinement to 
extend the spatial coherence in the ID geometry. The 
finite momentum modes are most probably excited due 
to dynamical nature of the acoustic lattice that interacts 
with the zero momentum polaritons, populated indirectly 
in the OPO setup, and scatters them to finite momenta. 
Therefore, for example, etching the whole planar struc- 
ture to form a wire will remove the main source of de- 
phasing that limits currently the coherence in the spatial 
domain. Such structures were recently produced [12] and 
a coherent propagation over a long distance of tenth of 
micrometers in these structures was reported [l2^ . 

Acoustic lattice induced regime of a ID condensate, 
which was identified in Section 6, opens a new way to 
study distinct non-equilibirum properties of the polari- 
tons. A better insight can be obtain in more complex 
frameworks, such as [30- 32] , that can capture finer de- 
tails and can provide a further understanding of the non- 
equilibirum processes in microcavities. 
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